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ABSTRACT 

The launch of the Hinode satelhte has allowed unprecedented high resolution, 
stable images of solar quiescent prominences to be taken over extended periods 
of time. These new images led to the discovery of dark upflows that propagated 
from the base of prominences developing highly turbulent profiles. As yet, how 
these flows are driven is not fully understood. To study the physics behind this 
phenomena we use 3-D magnetohydrodynamic (MHD) simulations to investigate 
the nonlinear stability of the Kippenhahn-Shliiter prominence model to the mag- 
netic Rayleigh- Taylor instability. The model simulates the rise of a buoyant tube 
inside a quiescent prominence, where the upper boundary between the tube and 
prominence model is perturbed to excite the interchange of magnetic field lines. 
We found upflows of constant velocity (maximum found 6kms~^) and a maxi- 
mum plume width ^ 1500 km which propagate through a height of approximately 
6 Mm in the no guide field case. The case with the strong guide field (initially 
By = 2Bx) results in a large plume that rises though the prominence model at 
~ 5kms~^ with width ~ 900 km (resulting in width of 2400 km when viewed 
along the axis of the prominence) reaching a height of ~ 3.1 Mm. In both cases 
nonlinear processes were important for determining plume dynamics. 

Subject headings: Instabilities — Magnetohydrodynamics (MHD) — Methods: 
Numerical — Sun: Prominences 
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Introduction 



From supernova outflows f lHachisu et aljl992l) t o the Earth's ionosphere ( iTakahashi et al 



20091 ) the Rayleigh- Taylor instability ( lTaylorlll950l ). in both its hydrodynamical (HD) and 
magnetohydrodynamical (MHD) forms, is a fundamental process in astrophysics. In Solar 
physics, the magnetic Rayleigh- Taylor instability is know to be an important driver behind 
many observed phenomena. For example the flux emergence process, where magnetic flux 
rises from the solar interior to the surface and undergoes the magnetic Rayleigh- Taylor in- 
stability as it eme rges into the upper atmosphere producing observed filamentary structure 
dlsobe et al.lbooeh . 



The physical properties of this instability have been well studied. The growth rate (oj) of 
the magnetic Rayleigh- Taylor instability for a uniform magnetic field parallel to the interface 



Att{p+ + p.)kgj 

where B is the magnetic field strength and k\\ is the component perturbation in the direction 
of the magnetic field and A is the Atwood number, where A = (p+ — p_)/ (p+ + p_) with p_|_ 



( Chandrasekhar 


1961 


). In the 


spikes is common ( 


Daly 


1967) 



in the magnetic field across the contact discontinuity, finding that this suppres sed the small 
wavenumbers creating wider filamentary structures. Recent observations by iBerger at al. 



( 120101 ) of dark upflows that propagate from underdense bubbles through quiescent promi- 
nences appear to be the observational signature of the Rayleigh- Taylor instability in quiescent 
prominences. 



Quiescent prominences are large structu res of relatively cool (10000 K iTandberg-Hanssen 



19951 )). dense (~ 10^^ cm~^ iHirayamal Il986l ) plasma, that exist in quiet regions of the solar 
corona, predomin antly at high he liographic latitudes. Using a chara cterist i c gas pressure 
of 0.6 dyn cm'"^ (IHirayamal Il986l ) and magnetic field of 3 ~ 30 G (jLeroyl Il989l ). gives a 
plasma /3 ~ 0.01-1. Linear magnetohydrostatic modelling of a quiet region filament has 
shown plasma /3 < 1 and strong departures from force-free magnetic field (jPudik et al. 



20081 ). Globally , quiescent prominences are incredibly stable structures that often exist in 
the corona for weeks. In contrast to this global stability, locally quiescent prominences are 
highly dynamic phenomena. Observations of quiescent prominence s have shown downflows 
( lEngvoldlll98ll ). vortices of approximately 10^ km xlO^ km in size ( iLiggett and Zirinlll9841) 



and a bubble of size 2800 km forming a keyhole shape with a bright center ( Ide Toma et al. 
20081 ) with velocities of 10-30 km s~^. There are many reviews that give a full description of 
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the curre nt understanding; of the structure and dynamics of quiescent prom inences (see, for 
example, lTandberg-Hanssenlll995l : iLabrosse et al.ll20ld : iMackay et al.ll2010l ) 



Observations by the Solar Optical Telescope ( iKosugi at al.ll2007l ) on the Hinode satellite 
( iTsuneta et al.ll2007l ) have show n that on a srnall sca le qu iescent prominences are highly dy- 
namic and unstable phenomena. iBerger et al.l (120081 ) and iBerger at al.l ( l2010l ) reported dark 
plumes that propagated from large bubbles (approximately 10 Mm in size) that form at the 
base of some quiescent prominences. Plumes form at the bubble prominence boundary and 
then flow through a height of approximately 10 Mm before dispersing into the background 
prominence material (see Figured]). Observations imply th at the plumes and the cavities 
have a column density about 20% of the prominence density (IHeinzel et al.ll2008l ). The dark 
upflows maintained an almost constant velocity of approximately 20kms~^ throughout their 
rise phase. Often these plumes would se parate from t h e larg e scale bubble forming smaller 
bubbles inside the prominence material. IBerger et al.l (120111 ) presents observations of large 
scale prominence bubbles using the Atmospheric Imaging Assembly (AIA) on the Solar Dy- 
namics Observatory (SDO) that show the temperature of the material inside the bubble to 
be > 250, 000 /ST. 



Berger at al.l ( 120101 ) hypothesized that the observed upflows were caused by the Rayleigh- 
Taylor instability. The observed dynamics seem to be those of a mixed mode perturbation 
with an interchange mode where a component of k, the perturbation wave vector, is perpen- 
dicular to the magnetic field B and a n undular mode wher e a component of k is parallel to 
B, in the high Atwood number limit. iRyutova et al.l (120101 ) described how the theoretically 
predicted growth rates and behavior for the magn etic Rayleigh- Taylor instability well match 
the observations of quiescent prominence plumes. Ivan Ballegooijen fc Cranmerl (120101 ) spec- 
ulated that the occurrence of the Rayleigh- Taylor instability inside a quiescent prominence 
could lead to the formation of the thread like structure observed through creation of a tan- 
gled magnetic field. If these observations are created by the Rayleigh- Taylor instability, it 
has very important implications for our understanding of quiescent prominences as it implies 
that gravity will play a very important role in the evolution of the system. 



The model th at we use in this work is the Kippenhahn-Schliiter model (iKippenhahn fc Schliiter 



19571 : |Priestlll982l )). The Kippenhahn-Schliiter model for solar prominences uses the Lorentz 
force of a curved magnetic field to support plasma against gravity (see Figure [2]) where 
magnetic tension balances gravity and gas pressure balances magnetic pressure. This model 
describes the local structure of the prominence, without including a corona, and is uniform 
in the vertical direction. This model has been shown to be linearly stable to ideal MHD 
perturbations (IKippenhahn fc Schliiterl Il957l : lAnzer I Il969l ) . 



Though this model only describes the local prominence geometry, neglecting the global 
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structure, it can be used to proved an important first step toward modeling these upfiows in 
a prominence like magnetic geometry to understand how the structure of the magnetic field 
evolves through time, as this can lead to a greater understanding of the quiescent promi- 
nence/filament system. In this paper we present a study of how a nonlinear perturbation 
in the form of a low density tube placed inside the Kippenhahn-Schliiter model can be per- 
turbed to allow the interchange of magnetic field lines causing upfiows to form inside the 
prominence model. 



2. Numerical Method 

In this study, we use the 3D conservative ideal MHD equations. Constant gravitational 
acceleration is assumed, but viscosity, diffusion, heat conduction and radiative cooling terms 
are neglected. Though the p artial ionisation of quiescent prominences may have some effect 



on the magnetic geometry (IHillier et al.ll2010l ). as the timescales of interest here are two 



orders of magnitude smaller than those associated with the effects of partial ionisation so 
they have been neglected in this paper. We assume the medium to be an ideal gas. The 
equations are non-dimentionalized using the sound speed (Cg = 13.2 km s~^), the pressure 
scale height (A = Cs/{'~fg) = RgT/{fig) = 6.1 x 10"^ cm), the density at the centre of the 
prominence {p{x = 0) = 10"^^ g cm~^) giving a characteristic timescale of r = A/C<i = 47 s. 
The value for 7 = 1.05 is taken, this is for simplicity as the internal energy e = p/{'-f — 1) 
excludes the possibility of using 7 = 1. Using /3 = 0.5, which gives an Alfven velocity of 
Va = Cs^/2/^(3 = 25.8 km s-^ 

The initial model is as follows: 

B,{x) = B^o (2) 
B,{x) = B,aanh(^^] (3) 



. i|.eo..- (i|^£) (4) 
where Bro is the value of the horizontal field at x = and B^oo is the value of the vertical 



field as X — )■ 00. This model is the Kippenhahn-Schliiter model as presented in iPriestI ( 1l982l ). 



As the Kippenhahn-Schliiter model is linearly stable to ideal MHD perturbations, a 
nonlinear perturbation is necessary. The perturbation considered here is a high temperature, 
low density tube placed in the center of the Kippenhahn-Schliiter model. The density of the 
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bubble at X = 2 = is 0.3p(0) (Atwood number A = 0.53) with width 2A and height 
8A. Figure [2] shows a visual representation of the initial conditions. The color contour 
represents the mass density, the lines represent the magnetic field lines. The initially the 
physical quantities are uniform in the y direction. To excite the interchange mode a velocity 
perturbation in Vy was imposed, where Vy was given as a sum of sinusoidal curves of different 
wavelength. The maximum amplitude of the perturbation {\vy\) is less than 0.01(7^. 

To reduce computational time, we assume a reflective symmetry boundary at x = 0. 
Due to the nature of the magnetic fleld at the top and bottom (z) boundary and at x = L^, 
the choice for boundary is very limited. A free boundary is assumed a.t x = Lx with a 
damping zone (damping time r = 4.4) for the hydromagnetic variables and Bz (to maintain 
the angle of the magnetic fleld at the boundary). This condition is important to maintain 
tension at the boundary to make sure the plasma can be supported there, but as this is a 
local simulation this will ultimately have an influence the results. For the top and bottom 
boundary, a periodic boundary is assumed and a reflective symmetry boundary is imposed 
at y = 0, Ly. 

The scheme used is a two step Lax-Wendroff scheme based on the scheme presented 
in lUgail (120081 ). using the artiflcial viscosity and smoothing also presented in this paper. 



The grid size is uniform in the y-direction, and in the x-z plane we take a grid of 75 x 400 
grid points , where the total area of the calculation domain is 3.5A x 85A. A flne mesh is 
assigned to an area of 40 x 320 grid points, actual size 1.2A x 30A, around the upper contact 
discontinuity allow the plumes to be resolved. In the y-direction a total of 150 grid points 
were used with dy = 0.05. 



3. Evolution of the Upfiows 

Figure [3] shows the evolution of the upflows for the simulation presented in this paper 
(where B.j.o = B^oo so (3 = Pgas/Pmag = 0.5). Upflows of size 2A-3A in width with velocities 
~ 0.39Cs can be clearly seen in panel D of the flgure. 

The evolution of the upflows can be understood in the following way. First the buoyant 
tube rises inside the prominence as described in the previous section. Then as the rise of 
the tube halts the interchange of magnetic fleld lines is excited by the small perturbation 
given to the system at the start of the simulation. The upflows excited have a wavelength of 
< 0.5A. This most unstable wavelength is decided by the numerical viscosity of the scheme 
and the terms added to the scheme for stability. 

As the upflows grow they interact with each other to create larger plumes. This interac- 
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tion is driven by the slight difference in plume size created by the random perturbation giving 
vortexes of different strength and orientation on either side of each plume. The interaction 
of the plumes drives stronger vortex motion which creates greater interaction between the 
plumes. The result of this interaction is the formation of the large plumes that dominate 
the system. This is known as the inverse ca scade process and is a common feat ure of the 



Rayleigh- Taylor instability (see, for example, lYoungsl (Il984j ) or llsobe et al.l ( 120061 )). 



Panels C & D of figure [3] show the rising bubble and falling spike formation associated 
with the high Atwood number Rayleigh- Taylor instability. Though the Atwood number in 
this case is only 0.53, due to the significant density difference, these dynamics are observed. It 
can be expected that they will become more pronounced as the Atwood number is increased. 
It should also be noted that the density difference and magnetic field suppress the role of 
the Kevin- Helmhotz instability. 

After the initial acceleration, the individual upfiows all achieve a relatively constant 
upflow velocity. This implies that at the top of the plume a force balance is maintained. 
This happens as the initial reduction in tension force from the upfiows is balanced by the 
accumulation of magnetic field at the head of the plume. 

The 3D structure of the magnetic field evolution caused by the upfiows and downflows is 
displayed in Figure IH The figure shows the density isosurface at p = 0.85 and the magnetic 
field lines at t = 15.3 & 43.6. Compared to the initial conditions, it is clear that the initial 
rise of the cavity reduces the tension in the magnetic field around the contact discontinuity. 
The figure shows that the upfiows form tubes inside the prominence and that the changes in 
the magnetic field distribution are small. The rise and fall of the field lines is initiated. The 
field lines move by gliding past each other in an interchange process. As the system evolves, 
the curvature of the magnetic field lines remains approximately constant after the distortion 
caused by the initial rise phase. 

From this the upfiows created inside the Kippenhahn-Schliiter model can advect the 
magnetic field further than the difference between the height of a magnetic field line at 
X = and x = xmax- This process can be seen as analogous to the interchange mode of 
the magnetic Rayleigh- Taylor instability where the position of magnetic field lines swapped 
without giving any change to the direction of the magnetic field. 

Figure [5] shows the evolution of the instability for the same initial conditions as above, 
but with an initial guide field of By = 23^0 giving a plasma (3 = 1/6. Due to problems with 
numerical reconnection in current sheets inside that are created in the underdense tube by 
shear to the guide field, initially the system is allowed to rise in 2D [x- z plane) without a 
guide field, and once the area around the contact discontinuity has reached an approximate 
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equilibrium this is taken as the initial conditions for the simualtion. The guide field is then 
added and perturbations to this system are solved. 

The upper 3 panels of figureOshow the 2D slice in the a; = plane. The initial formation 
of three plumes of ~ 1.3 Mm in size which through nonlinear coupling produce fiows of larger 
scale. The instability takes longer to grow than previously and the turbulent fiow of the first 
case is not present. The 3D images show that the large plumes seen in the 2D images are a 
result of cutting the filamentary structure at an angle, with the filaments aligned with the 
magnetic field. The rise velocity reaches 0.37 for a plume of width ~ 2.4 Mm in the y-z 
plane (~ 0.9 Mm in the direction perpendicular to the magentic field) with height ~ 3.1 Mm. 

4. Discussion 

In this paper, results of the nonlinear evolution of the magnetic Rayleigh- Taylor insta- 
bility in the Kippenhahn-Schliiter model are presented. The results show that the nonlinear 
effects are very important in determining the dynamics of the upflows. We found that nonlin- 
ear mode coupling was important for inducing an inverse cascade forming large upflows from 
smaller upflows. The dynamic pressure of the plumes was shown to be able to deform the 
magnetic fleld to some extent, but more importantly it allowed the interchange of magnetic 
field lines. Therefore, even though the magnetic field did not undergo severe deformation, 
upfiows were be excited. 

The model presented here is a local model of a prominence, where the global prominence 
structure can only be inferred through the boundary conditions. As this is a local simulation, 
it must be stressed that the boundary conditions play a strong role in the determination of 
the system. Also, it has been assumed that the magnetic field penetrating the hot tube is 
initially has the same distribution as the cool dense material. Even with these assumptions, 
a qualitative comparison between the prominence observations and the plume structure can 
be performed. For the comparison with these simulations and observations, the prominences 
plumes presented in Figure [T] are used. For easy comparison, we will assume that in the 
observations the magnetic field is aligned along the line of sight. 

Initially, the bubble forms beneath the prominence and then becomes unstable in mul- 
tiple areas. The perturbations form small-scale upfiows, which drive the creation of larger 
scale upfiows of a few megameters in width. These upfiows develop highly turbulent profiles. 
These upfiows described and the upfiows found in the simulations presented are very similar 
in appearance. 

To provide a more quantitative comparison the velocity and characteristic sizes can 
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be used. The observed upflows were found to supersonic, sub Alfvenic flows witli almost 
constant velocity of approximately 20km s~^. The simulations develop flows of an almost 
constant velocity with average velocity of 5.1 kms~^. This is more than a factor of 3 smaller 
than the upflows observed. The observed plumes propagate through a height of approxi- 
mately 10 Mm and have a characteristic width of ~ 300 km— 2 Mm. The simulations pro- 
duced upflows that have an initial width of ~ 200km, but through nonlinear processes 
produced upflows of ~ 600 km— 1.8 Mm in width and in case 1 at the end of the simulation 
had propagated through a height of 6 Mm. 

The inclusion of a guide field does not change the dynamics drastically, where the main 
differences are that the instability takes longer to grow and creates much larger structures 
in the y- z plane. The smaller growth rate could be a result of the stronger field strength 
suppressing the perturbation. The increase in field strength is likely to be the cause for the 
visible reduction in the turbulent like profile, as in this case the magnetic pressure dominates 
the dynamic and gas pressure. It will be important to investigate a guide field that creates 
shear between the magnetic field of the large bubble and the prominence , as the size of the 



filaments is likely to be dependent on the shear of the magnetic field (jStone fc Gardiner 



20071 1. 



In this paper, a boundary that keeps the magnetic field pointing upward at the bound- 
ary is used. This boundary confines the prominence material, chosen because of the lack of 
downfiows of prominence material through the large bubbles formed at the base of promi- 
nences due to upwardly curved magnetic field. A fixed boundary was not used due to the 
clear motion of plasma (and with it magnetic field). Therefore the boundary used was seen as 
the most probable. To elucidate this further, a full study of the boundary will be presented 
in a separate paper. 

The paper presents the first efforts to simulate the formation of upflows in a solar 
prominence. The observation of these prominence upflows has presented many difficult 
questions, among which are 'how do turbulent flows develop in the low-beta environment 
of prominences?' and 'what is the origin of the upflow plumes and how do their dynamics 
depend on factors such as bubble density and prominence magnetic fleld strength?' The 
answers to these questions will have implications not only for prominence science, but for 
the study of complex flows in any low-beta magnetic plasm a environment, for ex ample the 



Equatorial Spread-F phenomenon in the Earth's ionosphere (ITakahashi et al.ll2009l ). To fully 
understand the dynamics of these prominence bubbles and the plumes generated therefrom 
will require a wide parameter survey using the model developed here. 



Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as 
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Fig. 1. — Observations showing the formation of dark plumes in a quiescent prominence 
observed at 8-Aug-2007 20:01 UT observed in the Ha hne. The axes are labelled in mega- 
meters. 




Fig. 2. — Contour plots of the initial density distribution for the standard model for A) the 
X — z plane at ?/ = (with magnetic field lines), B) the y — z plane at x = & C) the 3D 
visualisation. All physical quantities are initially constant in the y direction. The initial 
velocity perturbation is applied to the upper contact discontinuity in along the y direction. 
The lengths displayed are given in megameters. 
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Fig. 3. — Temporal evolution of upfiows for t = 719, 1452, 2049 & 2453 s (normalized units 
t = 15.3, 30.9, 43.6 & 52.2) taken in the y — z plane at x = 0. All lengths are given in 
megameters. 




Fig. 4. — 3d structure of bubbles created in prominence with seleceted magnetic field lines 
at times t = 719 & 2049 s (normalized units t = 15.3 & 43.6.) The online color version 
distinguishes the field lines that initially thread the hot tube (gray) and those that don't 
(rainbow) by color. 
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Fig. 5. — Panels A,B & C show the temporal evolution of upflows for t = 1023, 1646 & 
2435 s (normalized units t = 22.7, 36.6 & 52.1) taken in the y — z plane at x = for the 
case with a guide field of strength By = 2B,j.q. All lengths are given in megameters. Panels 
D & E show the 3D structure at times t = 1646 & 2435 s. 



